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Abstract. The results of a NLTE model atmosphere 
analysis of 27 hydrogen-rich central stars of old plan- 
etary nebulae (PN) are reported. These stars were se- 
lected from a previous paper in this series, where we 
gave classifications for a total of 38 central stars. Most 
of the analyzed central stars fill a previously reported gap 
in the hydrogen-rich evolutionary sequence. Our observa- 
tions imply the existence of two separated spectral evo- 
lutionary sequences for hydrogen-rich and -poor central 
stars/white dwarfs. This is in line with theoretical evolu- 
tionary calculations, which predict that most post-AGB 
stars reach the white dwarf domain with a thick hydrogen 
envelope of w IO^^'Mq. 

We determine stellar masses from the comparison with 
evolutionary tracks and derive a mass distribution for the 
hydrogen-rich central stars of old PNe. The peak mass and 
the general shape of the distribution is in agreement with 
recent determinations of the white dwarf mass distribu- 
tion. 

The properties of most analyzed stars are well ex- 
plained by standard post-AGB evolution. However, for 
eight stars of the sample other scenarios have to be in- 
voked. The properties of three of them are probably best 
explained by born again post-AGB evolution. Two of these 
are hybrid CSPN (hydrogen-rich PG1159 stars), but sur- 
prisingly the third star doesn't show any signs of chemi- 
cal enrichment in its atmosphere. The parameters of five 
stars are not in accordance with post-AGB evolution. We 
discuss alternative scenarios such as the stripping of the 
hydrogen-rich envelope by a companion during the first 
red giant phase or the formation of a common envelope 
with a possible merging of both components. Two stars 
(HDW4 and HaWeS) remain mysterious after all. They 
resemble ordinary hot DA white dwarfs, but due to very 
large evolutionary ages the presence of a PN cannot be 
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explained. We speculate that the nebulae may be shells 
produced by ancient nova outbursts. 

A wide spread of helium abundances is observed in 
the photospheres of central stars of old PNe. It is shown 
that a good correlation between helium abundances and 
luminosity is present. It is inferred that when the stars' 
luminosities fall below L « 300Lq depletion starts and 
the helium abundance steadily decreases with decreasing 
luminosity. The existence of this correlation is in qual- 
itative agreement with recent theoretical calculations of 
gravitational settling in the presence of a stellar wind. 

Key words: stars: abundances — planetary nebulae: gen- 
eral — white dwarfs — stars: fundamental parameters 



1. Introduction 

Central stars of planetary nebulae (PN) are the immedi- 
ate precursors of white dwarfs. Most stars enter the white 
dwarf cooling sequence through this evolutionary chan- 
nel (Drilling & Schonberner 1985| ). After the maximum 
effective temperature is reached (100,000 K and more) the 
nuclear burning ceases, the surrounding nebula disperses, 
and the central star stage ends. Thus the nuclei of old 
PNe mark the transition to the white dwarfs. During this 
stage the onset of gravitational settling, which causes the 
chemical purity of many white dwarf atmospheres, can be 
observed. 

Central stars of old PNe play a crucial role in our 
understanding of the formation of the two distinct se- 
quences of central stars of planetary nebulae (CSPN) 
and white dwarfs: the hydrogen-rich and the hydrogen- 
deficient (helium- and carbon-rich) one. For the hydrogen- 
deficient CSPN there exists a continuous sequence from 
the luminous central stars (Wolf-Rayet stars of spectral 
type [WC]) to the low luminosity DO white dwarfs via 
the PG 1159 stars, hot pre- white dwarfs with typical tem- 



peratures of 100000 K. Mendez (1991) presented a compi 



lation of spectral types for (mostly luminous) CSPN and 
found a fracti on of about 1/3 hydrogen-deficient objects. 
Liebert (1986) reported a ratio of H-poor to H-rich stars 
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of 1:(7±3) for the hot white dwarfs (spectral types DO 
and DA resp.) in the Palomar-Green survey (Green et 
al. |1986 ), but noted a lack of very hot DA white dwarfs 
(the counterparts of the PC 1159 stars). Holberg ( |1987D se- 
lected a sample of the apparently hottest DA/DAO white 
dwarfs from the PG survey and concluded that none of 
them appear to have a temperature in excess of 80000 K, 
with Teff typically in the range 60000 K to 70000 K. How- 
ever, detailed analyses of these stars were not presented. 



Fontaine & Wesemael (1987) discussed this seeming 



lack of very hot H-rich white dwarfs and proposed that 
all CSPN finally lose their hydrogen-rich surface layer 
and evolve through one channel: the PG1159/DO stars. 
After the stars have become white dwarfs gravitational 
settling causes small traces of hydrogen, hidden in the 
helium layer, to float up and form a very thin H-layer. 
This scenario is in con tradic t ion t o theoretical calcula- 
tions (e.g. Sc honbc rncr 1981 , 1983 ; Vassiliadis & Wood 
1994; Blocker 1995), which predict a "thick" H-layer of 
RilO"'*MQ, but could nicely explain many properties of 
white dwarfs known at that time (EUV/X-ray observa- 
tions, the DB gap, pulsational properties of ZZ Ceti stars). 
Due to improved observations and increased theoretical 
knowledge the interpretation has changed during the last 
years and it is now generally accepted that the vast ma- 
jority of DA white dwarfs possess a thick H-layer (cf. the 
recent review of Fontaine & Wesemael 1997). However, the 



interpretation of many observational results is not com- 
pletely unambiguous. Vennes et al. (1997), e.g., analyzed 
a sample of EUV selected DA white dwarfs and compared 
their results with cooling tracks calculated with thick and 
with very thin layers. They found better agreement with 
the latter tracks and interpreted this result as evidence of 
a very thin hydrogen layer of most DA white dwarfs (but 

1999| ). 



see discussion in Napiwotzki et al. 



The Fontaine & Wesemael (1987) scenario makes one 
prediction concerning the stars just entering the white 
dwarf cooling sequence: the vast majority of them should 
be hydrogen-poor PG1159 stars or DO white dwarfs. 
Although a small number of hydrogen-rich stars in this 
transition region don't contradict the very thin layer sce- 
nario (see discussion in Fontaine & Wesemael 1997| ), the 
hydrogen-poor ones should be the dominant species. The 
transition between CSPN and white dwarfs is represented 
by the central stars of old PNe. This prompted us to start 
a systematic survey of the nuclei of old PNe. We selected 
apparently very old, faint and extended nebulae whose 
central stars could be identified. Important sources were 
the article of Kwitter et al. ( 198^ ) and the compilation of 
nearby PNe by Ishida & Weinberger (1987). The central 
stars should be far evolved, just entering the white dwarf 
cooling sequence. First results on interesting stars were re- 



ported in Schonberner & Napiwotzki (1990; Paper I) and 
Napiwotzki & Schonberner (1991a; Paper II). A spectral 



turned out that 28 central stars have a hydrogen-rich sur- 
face composition and that only seven were hydrogen-poor. 
The ratio of hydrogen- rich and -poor objects in our sam- 
ple of central stars just entering the white dwarf sequence 
amounts to 4:1, in good agreement with the previous in- 
vestigations of Liebert ( 1986| ) and Mendez (1991) men- 
tioned above. Our result indicates that the ratio of H-rich 
and H-poor objects remains approximately constant dur- 
ing the evolution from the central star phase to the white 
dwarf cooling sequence, co ntrary to what is expected from 
the Fontaine & Wesemael ( 1987| ) scenario. Since this find- 
ing is of fundamental meaning for our understanding of 
many aspects of white dwarf evolution, we will present 
here the results of a model atmosphere analysis of our 
central stars. This will allow us to determine their precise 
evolutionary status. We will in this paper concentrate on 
the hydrogen-rich central stars, which are the key for our 
understanding of the pre-white dwarf evolution. Analyses 
of hydrogen-poor Wolf-Rayet CSPN were recently carried 
out by Leuenhagen et al. ( 1996| ) and Koesterke & Hamann 
(1997). A review on the properties of PG1159 stars is 



given by Dreizler et al. (1995a). Prior to our survey only 
three analyses of high gravity CSPN were published: A 7 
(Mendez et al. |1981| ), NGC 7293 (Mendez et al. |l988aD , 
EGB6 (Liebert et al. 1989| ). Our much enlarged sample 



classification of the complete sample of 38 CSPN is pro- 
vided in Napiwotzki & Schonberner ( |l995 ; Paper III). It 



now allows a more meaningful investigation of the stellar 
evolution in this region of the HR diagram. 

This paper is organized as follows. We start with a 
description of the observations and their analysis in Sect. 2 
and 3, respectively. The results of the model atmosphere 
analysis are presented and the evolutionary status of the 
CSPN is discussed in Sect. 4. The evolution of the helium 
abundance of hydrogen-rich central stars and white dwarfs 
is described and compared to theoretical expectations in 
Sect. 5. The article finishes with conclusions in Sect. 6. 
We will use the results presented in this article to discuss 
the question of the PNe distance scale in Paper V of this 



2. Observations and data reduction 

A description of the sample selection and observational 
details was already given in Paper III. Therefore we will 
provide here only a cursory overview and details on addi- 
tional observations. 

Spectroscopic observations were performed with the 
3.5 m telescope of the Calar Alto Observatory/Spain and 
the Cassegrain TWIN spectrograph. We performed three 
runs with low resolution set-ups (October 1989, Septem- 
ber 1990, and July 1991). Medium resolution follow-up 
observations were performed in July 1992, July 1995, and 
October 1995. The resolution of the low resolution spec- 
tra ranges from 5 to 8 A. Details on the first four runs are 
provided in Paper III. 

The July 1995 and the October 1995 runs used iden- 
tical set-ups. Grating T05 in the blue and T09 in the 
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red yielded a dispersion of 36 A/mm. Tek CCDs #11 and 
#12 with 1024 X 1024 pixels of size 22 fj.m were used. 
The resulting spectral resolution, as measured from the 
He and Ar lines of the comparison spectra, amounts to 
1.5. . . 1.8 A (FWHM). While we always exposed the region 
from 5600 A to 6800 A in the red, we needed two grating 
angle settings to cover the spectral range from 3700 A to 
5400 A in the blue chanel. However, we did not always 
get spectra for both ranges. The data reduction follows 
standard procedures and is identical to that described in 
Paper HI. Let us only remark that we tried to correct for 
the nebula (and night sky) emission lines. For PNe with 
moderate line strengths and smooth density distributions 
this worked well. In cases of PNc with high surface bright- 
ness and/or strong spatial inhomogenities contamination 
by nebular emission lines remained (either insufficient sub- 
traction or overcorrection). Examples are NGC6720 and 
HaWe 5. In these cases the cores of the stellar Balmer lines 
(Ha and H/3 most severely) and the He ii 4686 A line may 
be corrupted. 

We have added the CSPN PHL 932 and the hot DAO 



white dwarf HZ 34 from Bergeron et al. ( 1994 ) to our 
original sample (Paper III). On the other hand we didn't 
analyse the spectra of Sh 2-176 and IW2 from Paper III, 
because their quality is just good enough for a reliable 
classification but is too poor for a reasonable analysis. In 
the case of WeDe 1 (WDHS 1) we based our analysis on 
the spectrum taken by Liebert et al. (1994), because it's 
signal-to-noise ratio is higher than that of our spectrum. 

3. Model atmosphere calculations and fit 
procedure 

We calculated hydrogen and helium composed model at- 
mospheres with the NLTE code developed by Werner 
( |1986| ). Basic assumptions are those of static, plane- 
parallel atmospheres in hydrostatic and radiative equilib- 
rium. In contrast to the atmospheres commonly used to 
analyze DA white dwarfs, we relax the assumption of local 
thermal equilibrium (LTE) and solve the detailed statisti- 
cal equilibrium instead of the Saha-Boltzmann equations. 



As described in Werner (1986), the accelerated lambda it- 



eration (ALI) method is used to solve the set of non-linear 
equations. 

The analysis of white dwarfs is usually the domain 
of model atmospheres in LTE. It is assumed that devi- 
ations from LTE are kept small by the high densities in 
white dwarf atmospheres. Model atmospheres which drop 
the LTE assumption are somewhat loosely called non-LTE 
(NLTE) atmospheres. A detailed check of the validity of 
the LTE assumption was presented in Napiwotzki (1997). 
It could be shown that LTE is a very good approximation 
for the "classical" DA white dwarfs with temperatures be- 
low 40000 K. However, the situation is completely differ- 
ent for typical white dwarf central stars (Tcff = 100000 K, 
log 5 « 7.0). The higher temperature and lower gravity 



causes large deviations between both types of model spec- 
tra. Therefore it is necessary to use NLTE model atmo- 
spheres for the analysis of these stars. 

Levels and lines are included in NLTE up to n = 16 
for hydrogen, n = 32 for Hell, and n = 9 for Hei. Hei 
levels with n = 5 to 9 were merged into one singlet and 
one triplet level for every principal quantum number. The 
n = 5 levels were split with respect to their angular quan- 
tum number for the subsequent spectrum synthesis. Line 
blanketing by the Stark broadened lines is taken into ac- 
count consistently for hydrogen, Heii lines and the reso- 
nance series of Hei. The synthetic spectra are computed 
with the extended VCS broadening tables (Vidal et al. 
1970) provided by Schoning & Butler (1989; priv. comm.) 
and Lemkc (1997) for hydrogen and Heii. We used the 
He I line broadening tables computed by Barnard et al. 



( |1974D , Shamey ( |1964| ) and Barnard et al. ( |1969D . Profiles 
of the isolated He i lines 5016 A and 5876 A were computed 



according the prescription of Griem (1974) 



Pressure dissolution of the higher levels is described 
by the Hummer & Mihalas ( 1988| ) occupation probability 
formalism following the NLTE implementation by Hubeny 
et al. (1994). In contrast to our analysis of EUV selected 
DA white dwarfs (Napiwotzki et al. 1999| ) we did not in- 
crease the critical ionizing field /3crit adopted to calculate 
the occupation probability by a factor two as proposed 



by Bergeron (1993). The motivation of Bergeron was not 



a fiaw in the Hummer & Mihalas (1988) formalism, but 
a compensation for the inadequacy of the standard Stark 
broadening theory when line wings overlap. This improved 
the agreement with experimental measurements and the 
fits to DA white d warfs i n the Bergeron et al. ( 1992D sam- 
ple. Finley et al. (1997) reported also better fit quality 
for hotter EUV selected white dwarfs. However, the cen- 
tral stars analysed here have even higher temperatures, 
and their gravities are lower by typical one or two dex in 
comparison to the s tars discussed in Bergeron ( 1993 ) and 
Finley et al. ( 1997 ). It is by no means clear if the cor- 
rection factor is a valid approximation in the case of our 
CSPN. Thus we preferred to leave the original value of 
the critical ionizing field Pcrit unchanged. From our expe- 
rience with the EUV selected DA white dwarfs analysed 



by Napiwotzki et al. (1999) we estimate offsets of the or- 
der 0.1 dex in gravity and 2% for T^g when compared to 
an analysis, which applies this correction factor. 

Our NLTE model grids have to cover wide ranges in 
temperature (30000 K < T^s < 300000 K), gravity (4.75 < 
logg < 8.75), and hehum abundance (from pure hydrogen 
up to nuc/^-H = 3). Thus, to keep computational effort 
within reasonable limits, we didn't compute every model 
in this 3D-mesh, but only the needed ones. 

Due to the common lack of other temperature indica- 
tors both, effective temperature and gravity, must be de- 
termined from the Balmer lines. For DA white dwarfs this 
was successfully done by simultaneous line pro file fit ting of 
all available Balmer lines (e.g. Bergeron et al. 1992; Napi- 



4 



R. Napiwotzki: Spectroscopic investigation of old planetaries IV. 



X 
3 



> 



1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



Sh 2-216 



Ha 
HP 




H7 



KS 



83000 K 

44000 K 



-50 



AA/A 



50 



Fig. 1. Model spectra calculated for two different temper- 
atures are compared to the Balmer lines of LSV+46 21, 
the central star of Sh 2-216 



wotzki ct al. 1999). Generally the observed line profiles 
are well reproduced by model spectra with the optimum 
parameters. 

However, this method failed for nearly all hydrogen- 
rich central stars of old PNe (Napiwotzki 1992; Napiwotzki 
& Rauch 1994): no consistent fit to the Balmer lines (Ha 
to Hi5 were used) was possible. A strong trend is present: 
fitting of higher Balmer lines yields higher temperatures. 
This effect is demonstrated for the spectrum of the white 
dwarf central star Sh 2-216 in Fig. |l|. While a high T^s 
model (83000 K) fits the US line, only 44000 K are needed 
to fit the Ha line. This difference causes serious problems 
for the interpretation of the evolutionary status of white 
dwarf central stars. The previous investigation carried out 



by Mendez et al. (1988a) restricted the analysis to one 
Balmer (H7) line only and thus failed to recognize this 
discrepancy. 



Napiwotzki ( 1993a ) argued that the high temperature 
derived from the highest Balmer line available (H5 or He 
respectively) is closest to the real temperature of the white 
dwarf central stars. Important evidence came from the 
analysis of the sdO star BD-|-28°4211. Its spectral ap- 
pearance is very similar to central stars of old PNe. As 



can be seen from Fig. 3 in Napiwotzki (1993a) the Balmer 
line problem is of similar strength as it is in the central 
star of Sh 2-216 (Fig. 0). Since BD+28°4211 is relatively 
bright {V = 10.5) we were able to obtain a high-resolution 
spectrum and detected the 5876 A line of neutral helium 
(Fig. 4 of Napiwotzki 1993a). Hei lines are very sensi- 
tive temperature indicators in this parameter range. We 
derived a best fit with T^s — 82000 K in excellent agree- 
ment with the temperature determination from the He line 
(85000 K). 



Napiwotzki & Schonberner ( |1993| ) and Napiwotzki & 
Rauch ( |1994[ ) discussed possible reasons of the Balmer line 
problem: wind effects, magnetic fields, pressure ionization 
and line quenching, deficits in the line broadening theory, 
and modifications of the atmospheric structure due to line 
blanketing of heavy elements. The authors concluded that 
none of these effects is sufficient to explain the Balmer 



line problem. Kubat (1995) found a moderate sphericity 



effect on Ha and, to a lesser extent, on H/3, too, for pre- 
white dwarfs with log g « 6. These effects are considerably 
weaker for larger gravity and too small to account for the 
Balmer line problem. 

Finally, Werner ( |1996| ) presented the likely solution 
of this problem. He demonstrated that the inclusion of 
Stark broadening for C, N, and O lines can have a 
strong influence on the atmospheric structure of very hot 
hydrogen-rich stars. It is likely that the effect on the emer- 
gent spectrum is pronounced enough to solve the Balmer 
line problem. An application of his model atmospheres 
on BD+28°4211 gave good agreement with the observed 
spectrum, if the parameters of Napiwotzki ( 1993a[) were 
used. Werner's results confirm Napiwotzki's ( |l993a ) so- 
phisticated guess that the temperatures derived from the 
highest Balmer lines are the most reliable, close to the 
"real" temperature of the central stars. T he rea son for 
this behavior is evident in Fig. 1 of Werner (1996), which 



displays the temperature structure of model atmospheres 
computed with different treatment of line opacity. The 
inclusion of Stark broadened CNO lines causes a strong 
modification of the temperature structure in the atmo- 
spheric layers, where Ha and H/? are formed, while the 
change in the deeper layers, where US and He are formed 
is negligible. 

However, to apply this procedure the abundances of C, 
N, and O need to be known which is not the case for our 
sample. In addition the computation of a NLTE model 
grid accounting for the influence of C, N, and O in the 
above described manner requires very large amounts of 
computer time, hence we analysed the CSPN with our at- 
mospheres containing only H and He. Since the Balmer 
line fits result in the same gravity for every Balmer line 
(within the error limits) we used the following recipe: g 
was computed from the average of all Balmer lines and 
Toff was derived from the fit of HJ (and He if available) 
with g fixed at the average value. The temperatures de- 
rived in our early analyses (Napiwotzki & Schonberner 
1991b and partly in Napiwotzki 1993h| ) are often too low, 
because we could not use the US line at that time due to 
reduction problems. All our new analyses now include this 



line and the discrepant cases discussed in Pottasch (1996) 
with Zanstra temperatures much higher than our analysis 
results are now resolved. 

The line fits were performed with the least-square al- 



gorithm described in Bergeron et al. (1992). The observed 



and theoretical Balmer line profiles are normalized to a lin- 
ear continuum in a consistent manner. Wavelength shifts 
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are determined with a cross-correlation method and ap- 
plied consistently to each complete spectrum. The syn- 
thetic spectra are convolved to the observational resolu- 
tion with a Gaussian and interpolated to the actual pa- 
rameters with bicubic splines, and interpolated to the ob- 
served wavelength scale. 

The atmospheric parameters are then determined by 
minimizing the value by means of a Levenberg- 
Marquardt steepest descent algorithm (Press et al. 19921 ). 
Finally, an estimate of the internal errors can be derived 
fro m the covariance matrix. In contrast to Bergeron et 
al. ( 1992 ), we estimate the noise of the spectra (a) used 



for the X from the neighboring continuum of each 
hne. The S/N is adopted to be constant throughout the 
line. Line cores, which show some remaining contamina- 
tion from the nebula, were excluded from the fits. 

In our investig ation of EUV selected DA white dwarfs 
(Napiwotzki et al. 1999| ) we showed that the internal errors 
are usually much smaller than the real observational scat- 
ter as inferred from a comparisons of independent analysis. 



We used the data collection of Napiwotzki et al. (|1999|) to 
derive a more realistic estimate of observational errors. 
For this purpose we binned the data into 5000 K intervals, 
determined the scatter, and performed a simple linear fit 
with log Toff as independent variable. The results are 



CT(Toff) 



off 



-0.0995 -f 0.0273- log Tcff 



cr(log5) = -0.8673 + 0.2130 -log Toff 



(1) 
(2) 



A moderate extrapolation is necessary to apply these 
fits to our CSPN analyses. We took a conservative ap- 
proach and assumed that these a values correspond com- 
pletely to the external errors and added them quadrati- 
cally to the internal errors derived with the fit procedure. 

4. Results and evolutionary status of the central 
stars 

The results of our analysis are summarized in Table |l|. The 
fits to the Balmer lines and helium lines are displayed in 
Figs. U to 1^. Due to the Balmer line problem discussed 
above we derived usually different temperatures from the 
individual Balmer lines. The individual values are indi- 
cated in the plots. The spectra are grouped according to 
their spectral appearance. Within a plot the spectra are 
displayed in order of increasing galactic longitude. 

A Tas-g diagram with the results is given in Fig. ^j. For 
comparison evolutionary tracks of post-AGB stars com- 
puted by Blocker (1995), Schonberner (1982), Koester & 
Schonberner (198(:) labeled with the remnant mass are 



plotted as solid lines. The dashed low mass tracks are 
taken from the computations of Driebe et al. ( |1998D . They 
represent the evolution of stars whose hydrogen-rich enve- 
lope was stripped away by a companion during the (first) 
red giant branch phase (Kippcnhahn et al. 1967; Iben & 
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Fig. 3. Mass distribution of the central stars of old PN 
(solid histogram) compared to two mass distributions ob- 
tai ned fr om analyses of hot DA white dwarfs (Ber geron et 
al. |1992| , dotted histogram; Napiwotzki et al. |1999| , dashed 
histogram) 



Tutukov 1986). The remaining mass is too low for the 



ignition of helium burning and the star will eventually be- 
come a low mass white dwarf with a helium core. In the 
remainder of this paper we will term this scenario "post- 
RGB evolution" for the sake of simplicity. 

The dotted track describes the evolution of a 0.524Mo 



post-early AGB (post-EAGB) track from Blocker (|995|). 
This star suffers from two late helium flashes, which are re- 
sponsible for the extended loops, and eventually enters the 
white dwarf domain on a track very similar to that of the 
O.546M0 model. Helium burning stars with significantly 
lower masses will never enter the (early) AGB and thus 
never produce a PN (AGB manque evolution; cf. Dorman 



et al. 1993) 



The new analyses are plotted as filled symbols. Our 
new results are supplemented by data taken from the 
literature (open symbols) for hydrogen-rich CSPN and 
DA/DAO white dwarfs. Hydrogen-poor objects are rep- 
resented by PG1159 stars, DO white dwarfs, and 0(He) 
post-AGB stars. Stellar parameters and references are 
provided in Table A. 2 and Table |A.1| in the Appendix. 
The symbols of stars with a known PN are encircled. 
Although analyses of hydrog en-po or Wolf-Rayet CSPN 
exist (Koester ke & Hamann 1997), which indicate that 
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Table 1. Parameters of the hydrogen- rich central stars. Teff was derived from the fit of US (+He) as described in the 
text. The CSPN are denoted by their PNG designation and by their common names. Masses were derived from the 
comparison with the evolutionary tracks in Fig. 



PNG 


Name 


Tctt 


logs 


log TlHc /rlH 


M 


dNLTE 


R 








(K) 


(cms"""") 




(Mo) 


(pc) 


(pc) 


(lO'yrs) 


025.4-04.7 


IC 1295 


90100+6200 


6.66+0.30 


-1.31+0.12 


0.51+0.04 


715 


0.16 


7.6 


027.6+16.9 


DeHt2 


117000+6300 


5.64+0.22 


-0.79+0.12 




2382 


0.54 


27 


030.6+06.2 


Sh2-68 


95800+9300 


6.78+0.32 


-1.02+0.21 


0.55+0.03 


1054 


1.02 


200 


034.1-10.5 


HDWll 


68100+9400 


6.38+0.31 


-1.23+0.24 


0.39+0.03 


1176 


0.13 


6.6 


036.0+17.6 


A 43 


116900+5500 


5.51+0.22 


-0.14+0.22 




2649 


0.51 


25 


036.1-57.1 


NGC 7293 


103600+5500 


7.00+0.22 


-1.43+0.15 


0.57+0.02 


291 


0.69 


28 


047.0+42.4 


A 39 


117000+11000 


6.28+0.22 


-0.85+0.10 


0.57+0.02 


1931 


0.81 


22 


060.8-03.6 


NGC 6853 


108600+6800 


6.72+0.23 


-1.12+0.09 


0.56+0.01 


436 


0.42 


13.4 


063.1+13.9 


NGC 6720 


101200+4600 


6.88+0.26 


-1.14+0.09 
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these stars are counterparts of the lu minou s hydrogen- 
rich CSPN analyzed by Mendez et al. ( 1988a ) there is no 
simple way to place them into our Tcs-g diagram, because 
a g determination is not possible without further assump- 
tions. 

The central stars of old PNe fill nicely the re- 
ported "gap" of the hydrogen-rich sequence. A continuous 
hydrogen-rich sequence from the central star to the white 
dwarf region is revealed, and therefore evolutionary sce- 
narios claiming a hydrogen-poor stage of all (or most) pre- 
white dwarfs (Fontaine & Wesemael 1987 ) can be ruled 
out. The gap in the hydrogen-rich sequence can be traced 
back to selection effects. Some of the hottest hydrogen- 
rich white dwarfs in the Palomar- Green survey (Green et 
al. |1986D , the most important source of faint blue stars 
until recently, were originally misclassified as sdO or sdB 
stars (cf. Jordan et al. 1991, Liebert & Bergeron |l99|). 
It is much easier to pick out the peculiar hydrogen-poor 
PG 1159 stars and DO white dwarfs. In addition, previous 
analysis of very hot DA or DAO white dwarfs tended to 
produce too low temperatures. This is probably caused by 
the neglect of NLTE effects and the Balmer line problem 



in th ese stars (cf. the discussion in Liebert & Bergeron 
1995|) . 

We determined masses of the central stars by inter- 
polating in the evolutionary tracks in Fig. [4 The values 
are provided in the sixth column of Table |l|. Error lim- 
its were estimated by propagating the errors of temper- 
ature and gravity. Additional errors may result from the 
adopted tracks. This is not a major problem for stars lying 
on the post-AGB tracks. However, the mass determination 
of stars, which cannot be explained by standard post-AGB 
evolution, might be wrong, if the star didn't evolve accord- 
ing to the adopted scenario but has a completely different 
history. The individual cases are discussed below. 

The resulting mass distribution (binned over 0.05Mq) 
is shown in Fig. |^. For comparison we show mass distri- 
butions derived for hot DA white dwarfs (Bergeron et al. 



1992; Napiwotzki et al. 1999 ) scaled down to the number of 



analyzed CSPN. The Bergeron et al. mass distribution was 



redetermined with the Blocker (1995) and Driebe et al 



(1998) evolutionary tracks. The Napiwotzki et al. (1999) 



analysis already applied these tracks. The Bergeron et al. 
sample contains mainly relatively cool DA white dwarfs in 
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Fig. 2. T'off-.g diagram with the resuhs of the new analyses plotted as filled symbols. These data are supplemented 
by results taken from literature (open symbols). Hydrogen- rich objects are marked by squares, hydrogen-poor by 
triangles. Central stars are encircled, a dashed circle is used for suspected PNe. Since the parameters from literature 
are sometimes rounded to identical values the positions of these stars are shifted by small random values. Post-AGB 
evolutionary tracks are drawn as solid lines, labeled with the mass of the remnant. Dashed lines indicate the post-RGB 
tracks discussed in the text. The O.524Af0 post-EAGB track is drawn as dotted line 



the temperature interval 15000 K < Toft < 30000 K, while 
the Napiwotzki et al. sample is EUV selected and cov- 
ers 25000 K < Toff < 55000 K. The CSPN mass distribu- 
tion peaks in the 0.55 . . . 0.60 interval, similar to the white 
dwarf distributions. Note, however, that the low mass tail 
is more pronounced than in the DA studies. Napiwotzki et 
al. described that EUV selected samples do select against 
low mass white dwarfs, which explain the deficiency in 



their study. It is likely that slowly evolving low mass stars 
are preferentially detected in extended, low density nebu- 
lae. However, we don't expect a strong selection against 
high mass CSPN in our sample, because the evolution 
of massive central stars slows down and becomes com- 
parable to that of low mass stars, when low luminosities 
are reached. Blocker (1995) claims that high mass CSPN 
evolve even slower than O.GMq remnants in this phase. 
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However, this depends on the treatment of AGB mass loss 
(of. discussions in Blocker |1995| and Vassiliadis & Wood 



1994) 



Since the mean might be heavily biased by few objects 
with extreme masses, it is not very useful for the com- 
parison of white dwarf mass distributions. Finley et al. 
( |l997|) proposed a Gaussian fit of the mass peak as a ro- 
bust estimate of the distribution peak. Due to small num- 
ber statistics and/or the high fraction of low mass stars 
the width of the best fitting Gauss curve is relatively high: 
O.14M0. This is twice the value typically found in the 
DA investigations and casts some doubt on the fit qual- 
ity. However, we redid the fit with the width held fixed at 
O.O7M0 and computed a peak mass only O.OOIMq higher. 
This convinced us that the derived peak mass of O.SSMq 
is reliable. The corresponding value for the Bergeron et 



al. (1992) sample (O.SGMq) and for the Napiwotzki et al. 
(1999) sample (0.59Mo) are in agreement. Stasihska et al. 
(1997) determined the mass distribution of CSPN from a 



distant independent approach using the nebular H/3 flux, 
the angular radius, the expansion velocity, and the stellar 
V magnitude. This method is based on a simple model of 



the PN evolution (Gorny et al. 1997). Stasihska et al. find 



a narrow mass distribution with «80% of all objects in the 
range 0.55 . . . O.65M0. The median of the sample amounts 
to O.6OM0. However, this value is slightly dependent of the 
adopted nebula mass in their PN model. Although this is 
a completely different approach it is in good agreement 
with our results for the central stars of old PNe. However, 
Stasihska et al. find no central stars with masses below 
O.55M0. This might partly be explained by their selec- 
tion criteria, which prefer brighter, younger PNe. 

We have used the values of effective temperature, grav- 
ity, and mass given in Table and the apparent magni- 
tudes given in Paper III to calculate the distances dNLTE 
of the CSPN. The results are provided in Table |l| together 
with the radius R and the kinematical age ikin of the PNe 
computed from the nebular data compiled in Paper III. 
See Paper V for more details. 

4.1. The D AO stars 

Most central stars of old PNe cluster around Tcs ~ 
100000 K and \ogg « 7.0. These are the DAO central stars 
displayed in Figs. ^ and ||. Among this group the central 
stars of weh known PNe hke NGC 6720 (the Ring nebula), 
NGC7293 (the Helix nebula), and NGC 6853 (the Dumb- 
bell nebula) are found. Another member is the central star 
of Sh 2-216, the closest known PN. According to the defi- 
nition of the spectral type DAO the spectra are dominated 
by the hydrogen Balmer lines, but the He II 4686 A line and 
sometimes Hen 4542 A and 5412 A are present. Generally 
these DAO central stars compare well with low mass post- 
AGB t racks . That's in contrast with the results Bergeron 
et al. (1994) for field DAO white dwarfs. These stars are 
cooler and are probably descendents of the extreme hor- 
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Fig. 6. Our reanalysis of the fie ld DAO white dwarf HZ 34 
from the Bergeron et al. (1994) sample 



izontal branch (EHB), the so-called AGB-manque stars, 
which stay hot during their post-HB evolution and never 
enter the AGB. Thus we conclude that DAO stars can re- 
sult from both evolutionary paths. In Sect. |^ we will show 
that, according to the observational evidence, the DAO 
phase is simply a transition stage passed by stars entering 
the DA sequence. 

Note that our re-analysis of HZ 34 yields a temper- 
ature 10000 K higher than the value given by Bergeron 
et al. (1994). The new parameters are consistent with a 
post-AGB nature. Actually this star has parameters very 
similar to other DAO central stars with a nebula. There- 
fore it is tempting to speculate about an as yet undetected 
PN around this star. However, a deep survey for nebula 



emission performed by Tweedy & Kwitter (1994) ended 
up with a negative result for HZ 34. Another star with 
parameters similar to the DAO central stars is the well- 
known sdO star BD+28°4211. The analysis presented in 
Napiwotzki ( |l993a| ) yielded T^s = 82000 K, logg = 6.2, 
and nue/nn = 0.1. These parameters place BD-|-28°4211 
in a region of the HR diagram, which is consistent with a 
star either close to the lower mass limit of post-AGB evo- 
lution or with post-RGB evolution. Interestingly, Zanin 



Weinberger (1997) detected recently a very large, faint 



emission region with an angular diameter of 5° centered 
on BD-f28°4211, which may possibly be the PN of this 
star. However, more observational data is needed before 
one can decide whether this is a chance alignment or really 
the PN of BD-f28°4211. 

The coolest DAO white dwarf found in our sample is 
GD561, the central star of Sh 2-174. It has an effective 
temperature of only 69000 K and cannot be explained by 
standard post-AGB evolution (cf. Fig. The status of 
GD 561 was already discussed in Tweedy & Napiwotzki 
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Fig. 4. Spectral fits of DAO central stars. The stars are ordered by increasing galactic longitude 
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Fig. 5. Spectral fits of DAO central stars (cont.) 



( |l994|) . The existence of a PN provides us with const raints 
on possible evolutionary scenarios. Bergeron et al. ( 1994 ) 
proposed an origin from the extreme horizontal branch 
(EHB). The parameters of GD 561 can be reproduced by 
post-EHB evolution (e.g. Dorman et al. 1993). However 
the PN remains unexplained. Since the horizontal branch 
evolution lasts wlO* years, any nebula produced on the 
RGB would be dispersed in the meantime. The mass of 
the remaining hydrogen-layer (< O.OIA/©) of EHB stars 
is much too low to produce a relevant nebula at the end 
of the horizontal branch stage. Obviously this region of 
the HR diagram is populated by stars which are not of a 
post-EHB origin as well. 



GD 561 is most likely the outcome of close binary evo- 
lution. During the giant stage (on the AGB or RGB) of 
the central star precursor, it can fill its Roche lobe or 



even enclose the companion in its envelope (common en- 
velope evolution) . Especially the latter causes heavy mass 
loss and a modification of the stellar evolution. The re- 
sult is a hot star inside an expanding shell, which can 
mimic a normal PN (s ee e.g . the review of Livio 1993 ). 



Tweedy & Napiwotzki ( 1994 ) showed that the parameters 
of GD 561 can be explained by the track of a O.296Af0 star, 
which lost its envelope during a common envelope event 



on the RGB (Iben & Tutukov 1986). However, the evolu- 
tionary age (4 • 10^ yrs) predicted by these calculations is 
uncomfortably large compared with the kinematical age of 
Sh 2-174 (40000 yrs; T able |l|). The new post-RGB tracks 
of Driebe et al. (1998) provide a much better agreement: 
the O.414M0 track reaches the position of GD 561 within 
100000 years after the RGB is left. Note that the PN age 
estimate is rather uncertain because the expansion veloc- 
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ity has not been measured. A photometric search for cool 
companions in the infrared J, H, K bands (Napiwotzki 



1995) found a cool companion of GD 561 separated by 



4". Furthermore, the preliminary analysis revealed that 
the central star itself has an infrared excess of 0™5 in K 
indicating binarity. The system might in fact be a triple 
one. 



4-2. The peculiar central stars K2-2, HDWll, and 
PHL 932 

It appears that the CSPN K2-2 and HDWll (Fig. 
are lower gravity counterparts of GD 561/Sh 2-174. Their 
positions in the HR diagram (Fig. ^ fit into the post- 
RGB scenario. Judged from the comparison with the 
Driebe et al. (1998) tracks their masses are slightly lower 
(O.39M0 and O.SSMq for HDWll and K2-2, respec- 
tively) and their theoretical post-RGB ages lie in the 
range 100000 . . . 200000 yrs. While the kinematical age of 
K2-2 (65000 yrs) is consistent with that prediction, the 
age of HDWll appears to be too low (6600 years). Note 
that, again, no expansion velocity has been measured for 
HDWll and thus the age was calculated by using the 
canonical value of 20km/s. If the expansion is slower the 
age of HDW 11 might be higher. 



While the age determination of HDWll might be in 
marginal agreement with the post-RGB scenario, this is 
certainly not true for PHL 932. From Fig. ^ we read of a 
mass of O.28M0. The age corresponding to the position of 
PHL 932 is 3-10^ yrs, the PN age amounts to 7700 yrs only. 
Although both age determinations are subject to severe 
uncertainties the discrepancy is so large that this scenario 
can be ruled out. PHL 932 might be the result of a common 
envelope event of a star with a degenerate CO core as 
calculated by Iben & Tutukov ( 198E| ; see Iben & Livio 



1983 



for a review on common envelope evolution). Mcndez et 
al. (1988b) speculated that both stars of the binary might 
have merged during the common envelope stage. Note that 
our mass estimate is only valid if PHL 932 is a post-RGB 
star. 



4-3. The hot, high gravity central stars of EGB 1 and 
WeDel 

Two high gravity CSPN are distinctly hotter and more 
massive than typical DAG central stars: EGB 1 (0.65AfQ) 
and WeDe 1 (O.68M0). Their spectra don't show the Hen 
4686 A line indicating that their spectral type is DA. How- 
ever, due to the high temperature and gravity of these 
objects the upper limits, which can be derived from our 
optical spectra are not very stringent: nuc/nn < 0.02 for 
both stars. Thus the absence of helium lines can already 
be explained by a modest depletion. 



XI 
3 



> 

■i-H 




' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' 

HDW 11 

logg = e.3a log He = -1.23 



PHL 832 
lose = 5-93 log He 



35000 K 
3S0OOK 




-1.53 
Ha 

HP 



Hell 4606 
5876 



Hel 5016 



Hel 4471 



I I I T I 



I I I I I I I I I T I I I I I I I I I 

logg = 6.09 log He = -1.55 




Hel 5B76 



. Hel 4471 

I 10% 

. ... I .... I .... I ... . 



-50 



50 



AA/A 



Fig. 7. Spectral fits of the peculiar central stars HDW 11, 
PHL 932, and K 2-2 



4-. 4- DA central stars 

Three cooler DA white dwarfs are present in our sample, 
too (Fig. DeHt5 has similar parameters as the DAG 
white dwarf GD 561/Sh 2-174 although it is slightly hotter 
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(Toff = 76500K vs. 69100K). This places DeHt 5 closer to 
the post-AGB tracks (Fig. ||), but a non post-AGB na- 
ture of DeHt 5 is more likely. Several scenarios are possi- 
ble (see discussion of GD 561 above), but DeHt 5 compares 
well with the post-RGB evolutionary tracks of Driebe et 



al. ( |199^ ). The kinematical age of the PN (129000 years; 
Table |1|) is in good agreement with the post-RGB age 
(«100000years) estimated from the OAUMq track. 

The DA white dwarfs HDW4 and HaWe5 are more 
mysterious. Their spectral appearance (Fig. ^) and stel- 
lar parameters (Tea = 47300 K, logg — 7.93 and T^s — 
38100 K, logg — 7.58) resemble ordinary hot DA white 
dwarfs (cf. Fig. ^). However, the existence of PNe around 
these stars is very difficult to understand. The theoreti- 
cal post-AGB ages of 47000 K and 38000 K white dwarfs 
amount to « 3 • 10^ years, and « 6 • 10^ years, respectively. 
That is much older than any PNe found in our sample 
(cf. Table I) and one would not expect a PN to survive 



that long. Moreover, the nebulae of HDW4 and HaWeS 
do not appear to be extraordinarily old. Their kinematical 
ages amount to only 2000 and 3000 years and their surface 
brightnesses are not atypically low. Actually, the nebula 
of HaWe5 is bright enough to contaminate the cores of 
the Ha and H/3 lines in the central star spectrum (Fig. ^. 

We are also not aware of any binary evolution mod- 
els which can explain the properties of these two central 
stars and their nebulae, because no scenario is able to 
produce a remnant which reaches the positions of HDW 4 
and HaWe 5 within a reasonable time after departure ei- 



ther from the AGB or RGB (see e.g. Iben k Tutukov 1985 



Iben 1986; Iben & Tutukov 1993). The fundamental reason 
is that the time scales in the white dwarf domain are de- 
termined by the gravitational- internal energy of the stellar 



remnant (Blocker & Schonberner 1990). If the stellar evo- 
lution is abbreviated by a companion stripping off matter 
from the primary, the structure of the resulting remnant is 
less compact than that of a normal post-AGB star of the 
same mass. That means the gravitational-internal energy 
content is larger and cooling needs even more time. 

Since it seems impossible to explain the existence of 
a PN around these stars, we investigated the hypothesis 
that the nebula HDW 4 and HaWe 5 only mimic the ap- 
pearance of a PN. It is well known that white dwarfs in 
cataclysmic binaries can eject shells during an outburst 
caused by a thermonuclear runaway after they have ac- 
creted a certain amount of hydrogen-rich material from 
the companion, usually a late type main sequence star. 
The mass of the ejected material is of the order 1O~'*M0 
(cf. Table 5.7 of Warner 1995| ), much smaller than the 



mass of a typical PNe. Sometime after this nova event the 
expanding shell can become visible as spatially resolved 
object with a morphology not muc h diffe rent from plan- 
etary nebulae (see e.g. Slavin et al. 1995; Tweedy 1995). 
We inverted the Shklovsky formula (see e.g. Eq. 1 in Pa- 
per HI) and estimated shell masses of HDW 4 and HaWe 5 
of 10^^ Mq and 2 • 10~^Mo, respectively, consistent with 
a nova interpretation. 

However, nova shells usually have expansion velocities 



of 500km/s or higher (see Slavin et al. 1995 and references 



therein). We have obtained a spectrum of the Ha region of 
HDW 4 with a spectral resolution of 1.63±0.02 A (FWHM; 
measured from Gaussian fits of comparison lamp lines). 
We applied the same fit procedure to the nebular Ha line 
and the Nil 6548/84 A lines and got an average width of 
1.85 ± 0.06 A. If we ignore thermal broadening this pro- 
vides us with an upper limit on the expansion velocity of 
2?^exp < 47km/s. That corresponds to typical PN expan- 
sion velocities, and is much lower than the usual velocities 
observed in nova shells, even if one considers projection 
effects. 

Until now we did not detect any binary companion. 
The low resolution spectra of HDW 4 and HaWe 5 end at 
6800 A and don't display any signs of a cool companion. 
Ongoing mass accretion from a companion must be very 
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low if any, because otherwise the atmospheres of HDW 4 
and HaWe 5 would be contaminated by helium and heavy 
elements and would not show the pure hydrogen spec- 
trum we observe. A probable counterpart of our myste- 
rious stars might be the very old nova CKVul which ex- 



ploded in 1670. Shara et al. (1985) recovered it as a very 
faint star inside a faint nebulosity on deep CCD images. 
From the size of the nebulosity Shara et al. estimated an 
expansion velocity of 59km/s, lower than for any other 
known nova shell. However, a direct measurement of line 



widths is still lacking. Naylor et al. ( 1992 ) took the nova 
nature of CK Vul into question because the proposed cen- 
tral star doesn't show Ha in emission as it is usually ob- 
served even in ol d nova systems (Ringwald et al. 1996| ). 
Shara et al. ( 1985| ) concluded that any mass transfer in the 
CKVul system must virtually have stopped. An infrared 
investigation gave evidence for cool dust in the system 
(Harrison 1996 ). Although observational evidence is spu- 
rious it may well be that HDW 4, HaWe 5, and CK Vul be- 
long to the same class of objects. Theoretical calculations 
(Shara et al. |1993| ; Prialnik & Kovetz |1995[ ) predict that 
nova outbursts of white dwarfs with "atypical low masses" 
of 0.60 or O.65M0 (roughly the masses of HDW 4 and 
HaWe 5) should be very slow and eject high shell masses 
compared to nova with more "typical" white dwarf masses 
(O.9IM0; Webbink[l990|). This might explain the observa- 



tional properties of HDW 4 and HaWe 5 discussed above. 
However, no definite conclusions are possible without fur- 
ther observational data. A deep IR search for companions 
in both CSPN will be of special importance. 



post-AGB tracks, where only one other CSPN of our sam- 
ple is found. All other stars are more evolved. The com- 
parison with the evolutionary tracks yields post-AGB ages 
smaller than 1000 yrs, which is short compared to the kine- 
matical ages of their PNe (25000 yrs and 11000 yrs for 
A 43 and NGC7094, respectively). However, if we con- 
sider the error bars and uncertainties of age determina- 
tions both ages can in principle brought into agreement. 
On the other hand their position in the HR diagram, dis- 
tinct from the other central stars of old PNe, argues for 
the reality of this age discrepancy. The absence of a PN 
around HS 2324-^3944 (Tes = 130000 K; logg = 6.2; Drei- 
zler et al. 1996) is another argument for the reality of the 
discrepancy. A 43 and NGC 7094 are good candidates for 
bo rn aga in central stars (Schonberner 1979 , 1983 ; Iben et 
al. 1983), which returned to the AGB after a late thermal 



flash and repeat the post-AGB evolution, now as helium- 
burning stars. The result is a seemingly young central star 
in an old PNe. This explains the age discrepancy we found 
for the three hybrid stars discussed above. Fittingly, this 
scenario is often invoked to explain the chemical abun- 
dances of PG 1159 stars. However, we found another cen- 
tral star (DeHt 2) with very similar temperature and grav- 
ity and an old PNe (ikin = 27000 yrs). If the age discrep- 
ancy found for the hybrid central stars is explained by 
the born-again scenario it should apply to DeHt 2 as well, 
but no sign of a chemical enrichment of the atmosphere is 
present. This star might indicate that born again evolu- 
tion is not necessarily linked to strong chemical processing 
of the photosphere. 



4-5. Central stars of hybrid spectral type 

Our sample contains three hydrogen-rich PG1159 stars 
or hybrid CSPN (Sh2-68, A 43, NGC 7094). Their spec- 
tra show a Civ absorption trough in the region of the 
Hell 4686 A line as it is typical for the hydrogen-poor, 
carb on- a nd helium-rich PG1159 stars (see Dreizler et 



al. 1995a for a review). However, unlike in PG1159 stars 
the Balmer lines are visible, too, which indicates a much 
higher amount of hydrogen in their atmospheres. The first 
star of this type, Sh2-68, was found during our su rvey o f 
central stars of old PNe (Paper II). Dreizler et al. ( |l995aD 
analyzed A 43 and NGC 7094 taking also the lines of C, 
N, and O into account. They derived Tes = 1 10000 K and 
log 5 = 5.7 for both stars. The agreement with our analy- 
sis of A 43 is good, while we derive a significantly higher 
temperature for NGC 7094. The difference of the Balmer 
lines of both stars (Fig. indicates that the parameters 
of these stars cannot be equal. NGC 7094 is likely to be 
hotter and/or of lower gravity than A 43. HS 2324-^3944 
is the only star of this class without a PN (Dreizler et 
al. 1996| ). A search for a nebula around this star was not 
successful (Werner et al. 1997). 

The hybrid CSPN (hydrogen-rich PG 1159 stars) A 43 
and NGC 7094 are placed on the horizontal part of the 



4-6. Comments on individual objects 

PNG 025.4-04.7 (IC1295): The PN is relatively bright 
and produces severe contamination of the core of H/3 and 
the He 11 4686 A hue. Thus we excluded 4686 A from the fit 
and estimated the helium abundance from the He 11 lines 
at 5412 A and 4542 A. Only the wings of the H/3 line were 
fitted. 



PN G 036. 0+17.6 (A 43): The He 11 4686 A region is heav- 
ily contaminated by C iv lines and was excluded from the 
fit. 



PNG 036.1-57.1 (NGC 7293): The centra l star was pre- 
viously analysed by Mendez et al. ( |l988a[ ). Their result 
was Toff — 89000 K and logg = 6.9. The gravity is virtu- 
ally identical with our value logg = 7.00 ± 0.22, but their 
temperature is somewhat lower (Teg = 103600 ± 5500 K). 
This is easily explained when one takes into account that 
H7 was the only analysed hydrogen line. Our fit of H7 
results in Teff = 89000 K, in perfect agreement with the 
Mendez et al. result. 
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PNG 047.0+42.4 (A 39): The McCarthy et al. ( |1997| ) re- 
sult Tefi = 110000 logg = 6.3, and nnc/nn = 0.11 
agrees with our parameters. 



dependence of the fit result on the adopted helium abun- 
dance is much lower for our NLTE analysis than it is for 
LTE analyses (cf. Napiwotzki |997|). 



PNG 60.8-03.6 (NGC6720): The cores of H/3, H7, and 
Hell 4686 A are contaminated by the nebula lines and 
were excluded from the fit as well as the interstellar Call 



H and K absorption in He. McCarthy et al. (1997) fitted a 
high resolution spectrum obtained with the HIRES echelle 
spectrograph of the Keck I telescope and derived a some- 
what low temperature of the central star (80000 K). High 
resolution spectra are certainly useful if one has to deal 
with bright nebulae. On the other hand it is difficult to 
reduce the very broad Balmer lines of a \ogg — 7.0 star, 
which are distributed over several orders of the HIRES 
echelle spectrum, without deformation of the line profiles. 

PNG 066.7-28.2 (NGC7094): The Heii 4686A region 
is heavily contaminated by C IV lines and was excluded 
from the fit. 

PNG 125.9-470 (PHL932): The temperatures derived 
from the Balmer lines and the ionization equilibrium of 
helium differ by 3000 K. This is certainly caused by our 
neglect of metal line blanketing. The Balmer line profile 
differences of the model spectra and the observed spectra 
are in close agreement to the expectation from the theoret- 
ical comparison in Fig. 10 of Napiwotzki ( |1997 ). Mendez et 
al. (|l988bD derived Toft = 37000 ± 2000 K, log 5 = 5.5 ±0.2 
which is in marginal agreement with our results. Part of 
the difference might be explained by the use of different 
model atmospheres and fitting procedures. 



PNG 158.9+17.8 (PuWel): The McCarthy et al. ( |1997D 
derive T^s = 110000 if, \ogg = 7.0, and nne/nn = 
0.02. The temperature is slightly higher than our result 
(94000 K). See discussion of NGC6720. 



PNG 197.4- 06. 4 (WeDel): WeDe 1 was first announced 
as a very hot, high-gravity star by Napiwotzki & 
Schonberner ( 1991bD . A LTE analysis of the moderate 
quality spectrum available at that time resulted in lower 
limits on T^s ranging from 100000 K to 200000 K (depend- 
ing on the adopted helium abundance) and a lower gravity 
limit of log (7 = 8.5. A NLTE analysis using a better spec- 
trum and model atmospheres similar to that described 
in the present paper yielded in Toff = 217000 K and 



\ogg = 7.35 (Napiwotzki 1995). Independently, Liebert 
et al. (1994) performed a LTE analysis and derived Toff — 
99000 K . . . 165000 K, again depending on the adopted he- 
lium abundance. Since the quality of the Liebert et al. 
spectrum is superior to that of our Calar Alto spectrum, 
we used their spectrum for our analysis. The result is 
TefT = 141000 ± 32000 K and log<? = 7.53 ± 0.32. The 



PNG 215.5-30.8 (A 7): The McCarthy et al. ^997^ re- 
suh Toff = 100000 JC, logg = 6.6, and tiho/^h = 0.05 
agrees with our parameters. 

5. Helium abundances and the onset of diffusion 



A wide spread of helium abundances is observed in the 
atmos pheres of hydrogen-rich central stars. Venues et al. 
( 1988 ) have shown that helium sinks down due to gravita- 
tional settling in the atmospheres of hydrogen-rich white 
dwarfs. Diffusive equilibrium is reached fast and no vis- 
ible traces of helium should be present in these atmo- 
spheres. Therefore it was considered the most consistent 
physical picture of DAO white dwarfs that their atmo- 
spheres should consist of an extremely thin hydrogen layer 
on top of the helium envelope. The helium lines should 
form in the deeper, helium-rich regions. In the context 
of the Fontaine & Wesemael (1987) scenario of the spec- 
tral evolution of white dwarfs DAO stars are transition 
objects from the helium-rich DO to the hydrogen-rich 
DA sequence, which have just built up an extremely thin 
hydrogen-layer. The very thin layer hypothesis can be 
tested by detailed fits of the He 11 4686 A profiles. The for- 
mation depth of this line is quite different for stratified and 
homogeneous atmospheres: the line is much broader and 
shallower in the stratified case. Napiwotzki & Schonberner 
(|1993D demonstrated that the He 11 4686 A fine of Sh 2-216 



is well fitted by a homogeneous model atmosphere while 
a stratified atmosphere could be excluded. Figs. ^ and || 
show also good agreement between the profiles computed 
from our homogeneous models and those observed for He 11 
4686 A. Bergeron et al. (1994) extended this investigation 
to their DAO sample and found also that the atmospheres 
of most stars are not stratified. However, they presented 
one example, PC 1305—017, with a probably stratified at- 
mosphere, which might, indeed, be an object caught dur- 
ing the transitio n from DO t o DA spectral type. 

Napiwotzki ( |l993b| , |1995D and Bergeron et al. ( |1994D 
noted a correlation between helium abundance and grav- 



ity. However, Bergeron et al. (1994) did not consider it 
to be significant and stressed a stronger correlation with 
Toff. We search for correlations between the helium abun- 
dance and temperature, surface gravity and luminosity in 
Fig. n2l Wc combine our results with the DAO analyses of 



Bergeron et al. (1994) and the CSPN analysis of Mendez 
et al. ( |19884 |1992| ). AU stars of the Mendez et al. sam- 
ple (with the exception of the peculiar objects EGB 5 and 
PHL 932) are on the constant luminosity part of the cen- 
tral star tracks. Some stars exhibit helium enrichment by 
dredge up of nuclear processed material, but none shows 
a significant depletion below the solar value. If we ignore 
these luminous Mendez et al. stars, indeed a correlation of 
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helium abundance with temperature becomes visible (up- 
per panel of Fig. p^ ) , however, not a very good one. On the 
one hand the hottest stars WeDe 1 and EGB 1 have upper 
limits below the solar value and on the other hand several 
of the cooler stars are only mildly depleted. A better cor- 
relations is found between nuc/nu and the surface gravity 
(middle panel). However, in this case there are exceptions 
from the general trend present, too. A far better corre- 
lation is found between helium abundance and luminos- 
ity (lower panel). All stars, whether of post-AGB nature 
or not, can be explained by one trend line, although the 
scatter seems to be larger than expected from the observa- 
tional errors alone. For all stars with a luminosity higher 
than «3OOL0 the helium abundance is close to the "solar" 
value nuc/nn = 0.1 or above. For stars with lower lumi- 
nosity the helium abundance steadily decreases with de- 
creasing luminosity. Since effective temperature and sur- 
face gravity are related to luminosity it is not surprising 
that a correlation with temperature is found in a sample 
with roughly equal gravity (Bergeron et al. 1994 ), or a 
correlation with gravity in a sample with roughly equal 
temperature (Napiwotzki 1993b| ; [1995 ). Although intrin- 
sic scatter seems to be present in the correlation of he- 
lium abundance with luminosity its existence is remark- 
able. The observed abundances of heavy elements in hot 
white dwarfs usually display a more or less stochastical 
dependence on stellar parameters without any clear trend 



(cf. Chayer et al. |1995| ). 

Recent theoretical calculations performed by Unglaub 



& Bues (1998) show that the helium abundances observed 
in DAO stars can be explained if a reasonable mass loss 
rate is adopted. Mass loss counteracts the gravitational 
settling of helium and slows it down, if the strength of the 
stellar wind is high enough. The mass loss of hot stars is 
radiativ ely dr iven and a strong function of luminosity (e.g. 
Abbott 1982 ). Thus a correlation between luminosity and 
helium abundance is also expected on theoretical grounds. 
The calculations of Unglaub & Bues ( 199^ ) show that the 
static equilibrium abundances are not reached within typ- 
ical evolutionary time scales of hot white dwarfs. There- 
fore the stellar abundances are the result of a dynamical 
equilibrium between diffusion time scales and the evolu- 
tion of the star. This might explain why sometimes stars 
with similar parameters show different photospheric abun- 
dances. An example is given by the DA DeHt 5 and the 
DAO Sh 2-174. DeHt 5 has a slightly higher temperature 
and slightly lower gravity than Sh 2-174. Thus we would 
expect the helium abundance of DeHt 5 to be higher than 
in Sh 2-174, but the spectrum Sh 2-174 shows an easily 
detectable Heii 4686 A line (Fig. ^) while the Hue is be- 
low the detection limit in the spectrum of DeHt 5 (Fig. 
However, if the helium abundance is the result of time de- 
pendent diffusion it should depend on the previous evolu- 
tion of the stars. Thus stars with similar parameters but 
different evolutionary histories can show different abun- 
dance patterns. We conclude that today observational ev- 



idence and theoretical calculations of the evolution of sur- 
face helium abundance show semi-quantitative agreement. 
What is needed for an improved understanding are self- 
consistent diffusion calculations taking the stellar evolu- 
tion into account and an extension of the comparison to 
as many heavy elements as possible. 



6. Conclusions 

We presented a NLTE model atmosphere analysis of 27 
hydrogen-rich central stars of old PNe. The stars were 
taken from Paper HI in this series, where we provided 
classifications for a total of 38 CSPN. The analysis of the 
H-rich CSPN is hampered by the Balmer line problem: for 
most stars a consistent fit to all Balmer lines was not pos- 
sible. Generally the highest Balmer line yields the highest 
temperature. The likely solution of this problem was pre- 



sented by Werner (1996). He demonstrated that the inclu- 
sion of Stark broadening for C, N, and O lines can have 
a strong influence on the atmospheric structure of very 
hot hydrogen-rich stars. The lower Balmer lines are most 
strongly affected and their analysis with standard NLTE 
model grids yields unreliable temperatures. It is likely that 
the effect on the emergent spectrum is pronounced enough 
to solve the Balmer line problem. Since the computation 
of a NLTE model grid accounting for the influence of C, N, 
and O in the above described manner requires enormous 
amounts of computer time, we analyzed the CSPN with 
our atmospheres containing only H and He and derived 
Toff from the fit of US (or He if available). 

A previously reported gap in the hydrogen-rich evo- 
lutionary sequence is filled by our sample and, therefore, 
shown to result from a selection effect. Taking into ac- 
count the hot white dwarfs analyzed by Bergeron et al. 
( 1994 ) we have now a continuous sequence from the CSPN 
to the white dwarfs. The observational findings are well 
explained by a two channel scenario with a hydrogen- 
rich and a hydrogen-poor sequence. It is no longer nec- 
essary to inv oke an one channel scenario (Fontaine & We- 
semael 1987 ) to explain the spectral evolution of pre- white 
dwarfs. Note that our findings do not rule out that some 
white dwarfs change their spectral type from DO to DA 
(see di scussi ons in Fontaine & Wesemael 1997 and Napi- 
wotzki |l998l) . 

Stellar masses were computed from a comparison 
with evolutionary tracks and a mass distribution for the 
hydrogen-rich central stars of old PNe was derived. The 
distribution has a sharp peak centered at 0.55Mq This 
value and the general shape are similar to the results 
of previous analysis of hot DA white dwarfs (Bergeron 
et al. |1992| ; Finley et al. |1997| ; Napiwotzki et al. [1999D 



and of a recent distance independent investigation of PNe 



(Stasihska et al. 1997). A relatively large fraction of low 
mass central stars is present in our sample. This might 
be explained by our selection criteria which tend to prefer 
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the old and thin nebulae around slowly evolving post- RGB 
stars. 

Although the majority of analyzed stars is well ex- 
plained by standard post-AGB evolution, there exists a 
number of stars for which other scenarios have to be in- 
voked. The properties of three stars are probably best ex- 
plained by born again evolution. Two of these are hybrid 
CSPN (hydrogen-rich PC 1159 stars), but surprisingly the 
third star doesn't show any signs of chemical enrichment 
in its atmosphere. The third hybrid CSPN in our sample 
(Sh2-68) has parameters which place it among the group 
of normal DAO central stars. 

The fundamental parameters of five stars are not in 
accordance with a standard post-AGB evolution. They 
are likely the outcome of a close binary evolution. Several 
scenarios are possible. Four stars are reasonably well ex- 
plained by the evolution of stars whose hydrogen-envelope 
has been stripped off by a nearby companion during the 
first red giant branch phase. The remnants will eventually 
become low mass white dwarfs with a degenerate helium 
core. Although the post-RGB evolution is considerably 
slower than post-AGB evolution the existence of low den- 
sity nebulae around these objects is in accordance with 
theoretical time scales. The kinematical ages of the ob- 
served PNe are consistent with the evolutionary post-RGB 
age, except for PHL932, where this scenario is ruled out 
by a very large discrepancy between predicted post-RGB 
age and kinematical age of the nebula. Two stars (HDW 4 
and HaWeS) remain mysterious. They resemble ordinary 
hot DA white dwarfs with evolutionary ages of several mil- 
lion years. However, this large age is in clear contradiction 
with the presence of a nebula. Close binary evolutionary 
scenarios do not resolve this discrepancy, because the rem- 
nants need several million years after the possible ejection 
of a nebula to reach the position of HDW 4 and HaWe 5 in 
the HR-diagram. We speculate that the nebulae of these 
stars are produced by nova-like events. However, no direct 
observational evidence can be given in the moment. 

Na piwotzki & Schonberner ( |l993|) and Bergeron et al. 
( |l994| ) have shown that the atmospheres of DAO stars are 
chemically homogeneous, i.e. not layered. We have shown 
that a good correlation between helium abundance and lu- 
minosity is present. Depletion starts when the stellar lumi- 
nosity falls below L w SOOL©, and the helium abundance 
steadily decreases with decreasing luminosity. Although 
the scatter is probably larger than expected from the 
observational errors alone, the correlation is surprisingly 
tight considering the more or less stochastic abundance 
variations usually found in hot white dwarfs (Chayer et 
al. 1995). The existence of the correlation is in qualitative 
agreement with recent theoretical calculations of gravita- 
tional settling in the presence of a stellar wind (Unglaub 
& Bues |1998D. 
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Appendix A: A compilation of CSPN and white 
dwarf parameters from literature 

Parameter determinations of CSPN and hot white dwarfs and 
related objects are widely scattered through literature. Many 
sources are not easily accessible. For this reason we provide a 
table with parameters and references of the objects displayed in 
Fig. ^ This compilation may be useful for future investigations. 
The most obvious application is the selection of candidates for 
deep imaging to look for faint PNe. We tried to collect all data 
for central stars and pre-white dwarfs with temperatures higher 
than 30000 K from literature, which were analyzed with state- 
of-the-art techniques. Stars alrea dy a nalyzed in this paper (cf. 
Table |l|) are omitted from Table |A.2| . We used the same lower 
temperature limit for the white dwarfs, but did not aim for 
completeness for the cooler DA wh ite d warfs (Tctt < 60000 K). 
The results are provided in Table A. 2 The stars are grouped 
according to their spectral c lasse s and their temperature. Ref- 
erences are given in Table A.l. Objects in globular clusters 



were excluded from this compilation. 

While the DAO and DA analyses are usually based on LTE 
mo del a tmospheres (with the exception of the Napiwotzki et 
al. 1999 analysis) , N LTE analyses were performed for all other 
stars listed in Table [A.2[ Since deviations from LTE produce 
significant e ffects in the hottest D A an d DAO white dwarfs 
(Napiwotzki 1997] ; Napiwotzki et al. 1999), one should be aware 
of possible systematic effects. 
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Fig. 9. DA central stars. The synthetic spectra of the he- Fig. 10. Fits of the hybrid PG 1159 stars. The Heii 4686 A 
lium lines are computed for the 3a limit of the helium lines of A 43 and NGC7094 were not used for the fit. The 
abundance observed and theoretical spectra are only shown for com- 

parison purposes 
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Fig. 11. Fit of the high luminosity central star DeHt 2 
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Fig. 12. Photospheric helium abundance as function of ef- 
fective temperature, gravity, and luminosity (top to bot- 
tom). Our results (filled squares; hybrid CSPN: filled di- 
amonds) are supplemented with the DAO white dwarfs 
analysed by Bergeron et al. ( 1994 ; ope n squares) and the 
CSPN analysis of Mendez et al. ( 1988a ; stars). Upper lim- 
its are marked by arrows 



